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Abstract 

The thermodynamics of excited nuclear systems allows one to explore the second-order phase 
transition in a two-component quantum mixture. Temperatures and densities are derived from 
quantum fluctuations of fermions. The pressures are determined from the grand partition function 
of Fisher's model. Critical scaling of observables is found for systems which differ in neutron to 
proton concentrations thus constraining the equation of state of asymmetric nuclear matter. The 
derived critical exponent (3 = 0.35 =L 0.01, belongs to the liquid-gas universality class. The critical 
compressibility factor P c /p c T c increases with increasing neutron number. 
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Understanding the behavior of nuclear matter under extreme conditions of densityand 
temperature is one way to gain insight into two-component finite quantum systems jl[ 
The basic properties of nuclear matter (consisting of protons, Z, and neutrons, N) as a 
function of temperature, density, and proton fraction are described by the nuclear equa- 
tion of state (EOS). The dependence of the EOS on the proton fraction is essential not 
only in describing features of heavy-ion reactions but also in describing characteristics of 
astrophysical environments such as the dynamics of stellar collapse and supernovae as well 
as the formation and structure of neutron stars 1|. Characterization of the critical point 
of nuclear matter (temperature, density and pressure) as a function of the proton fraction 
provides important information about the nuclear EOS. 

In nuclear matter, the nucleon-nucleon interaction is similar to the van der Waals inter- 
action with short-range repulsion and long-range attraction. This analogy strongly suggests 
the occurrence of a nuclear phase transition similar to a liquid-gas phase transition [3 
Experimentally this can be probed with heavy-ion collisions around the Fermi energy. Un- 
like van der Waals fluids, nuclei are finite, two-component systems. Most of the divergences 
usually linked to a phase transition in macroscopic systems are washed out in these small 
systems 8|. Moreover, the additional degree of freedom which is related to proton and 
neutron concentrations makes the phase transition more complex 9[. The existence of a 
nuclear phase transition is currently the subject of much investigation via caloric curves 
0, [h]], critical exponents U|, negative heat capacities 12] and other observables Q, Q. 
However, because of the assumptions made in these studies, the system could not be located 
in the pressure-density-temperature space 151 ] . 

In this paper, we report the experimental quantum temperatures and densities of the 
fragmenting system calculated by the quantum fluctuation method for protons presented 



a 



in Refs. 



16 



17J. Since the protons represent the vapor part, the derived densities and 



temperatures refer to the vapor branch of the 'liquid-gas'-like instability region. In addition, 
pressures calculated through the grand partition function of Fisher's droplet model 18[ are 
also presented. It is shown that the present data contain a signature of a liquid-gas phase 
transition. Scaling of physical observables to their critical values displays universality, i.e. 
independence from the proton-neutron asymmetry of the source. In turn this implies that 
details of the second-order phase transition do not depend on the nature of the particles 
(either classical or quantum) nor on their interaction. All systems display universal scaling. 
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This is quite amazing in view of the fact that quantum methods are used to derive properties 
of nuclear systems, at variance with the approaches used for classical macroscopic fluids. 

The experiment was performed at the Texas A&M University Cyclotron Institute. Beams 
of 64 Zn, 64 Ni and 70 Zn at 35 MeV/nucleon were incident on targets of 64 Zn, 64 Ni and 70 Zn 



respectively 19|, |20|. The charged particles and free neutrons produced in the reactions were 



measured with the NIMROD-ISiS An detector array 21|. The granularity and excellent 
isotopic resolution provided by the array enabled the reconstruction of the quasi-projectile 
(QP) in mass, charge and excitation energy. The QP is the large excited primary fragment 
of the projectile following a non-central collision with the target which will subsequently 



undergo breakup. The Neutron Ball [22j provided event-by-event experimental information 
on the free neutrons emitted during a reaction. The QP source was selected by means of 



event-by-event cuts on the experimental data similar to those used in Refs. |23|, [24J with its 
mass (A = Z + N) restricted to be in the range 54 < A < 64. The excitation energy was 
deduced using the measured free neutron multiplicity, the charged particle kinetic energies, 
and the energy needed for the breakup (Q- value). 

Recently, the neutron-proton asymmetry of the source was identified as an additional 



order parameter in the nuclear phase transition 



25 



5M27|. Therefore, to investigate such a 
dependence the data were sorted into four different source asymmetry (m s = (N — Z)/A) 
bins ranging from 0.04 to 0.24 with bin width of 0.05. The mean m s values corresponding 
to these m s bins are 0.065, 0.115, 0.165 and 0.215. The four m s bins will be referred to with 
their mean values in the rest of the discussion. The effects of QP excitation energies on the 
thermodynamic quantities were investigated by further gating the data into nine bins, each 
1 MeV wide, in the range of 1-10 MeV/nucleon. 

The temperatures of the different selected QPs are calculated with the momentum 
quadrupole fluctuation method reported in Refs. 



17| . The momentum quadrupole is 



defined as 



p: 



Py where p x and p y are the transverse components of a given frag- 
ment's momentum. This quantity is zero on average in the center of mass of the equilibrated 
QP. The longitudinal component, p z , is excluded to minimize contributions from the colli- 
sion dynamics, which manifest in the beam direction. In this paper, we use protons which 
are fermions as our probe particle. Assuming a Fermi-Dirac distribution, the normalized 
transverse fluctuation of the fragment momentum quadrupole {(J xy /N) is connected to the 
temperature T by the relation 
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°f = Am 2 T 2 FQ C , (1) 

where N and m are the average multiplicity per event and the mass of a given fragment 
type, respectively. Fqc is the quantum correction factor which should converge to one at 
high T (classical limit) and is expressed as 

fTY in 

F QC = 0.2 (-J + 1 , (2) 

where £f = £/ (p/Po) 2 ^ 3 is the Fermi energy of the nuclear matter at density p. The quan- 
tities Pq and Ef denote the normal nuclear density and the corresponding Fermi energy, 
respectively. The values of po = 0-15 fm -3 and £/ =36 MeV are used in the calculations. 

A similar derivation is given for the normalized multiplicity fluctuation of fermions emit- 
ted from the QP, and it is shown to depend on T/ef. This quantity in turn is parametrized 
in terms of a 2 N / N which is given as 

T , 0.442 
— = -0.442 + 

e f (\ - & 

+0.345^-0.12^ . (3) 

Once the normalized fluctuations are experimentally determined for a given excitation en- 
ergy, using Eq. [2] and Eq. |3]the quantum correction factor Fqc can be obtained. From Eq. CEJ 
one can easily calculate the temperature T and, from the fact that £/ = £ / (p/ Po) 2 ^ 3 , the 
density p is obtained. More details can be found in Ms. HQ. 

Temperatures and densities of the QP are extracted through the momentum quadrupole 
and multiplicity fluctuations using protons as a probe particle. These are plotted as a 
function of excitation energy per nucleon in Fig. [TJ The protons represent the gas or the 
low-density region in the liquid-gas type phase transition. The extracted densities for the 
four source asymmetries (bottom panel of Fig. [T|) show a dependence on the value of m s . 
An overall ordering in the density with m s is observed for each excitation energy: the larger 
the asymmetry, the lower the density. A previous analysis has shown an ordering of the 
temperatures within a classical treatment 28|]. In the present treatment, the correlation 
between temperatures and densities indicates that the dependence on m s is manifest in the 
densities. However, the spacing between the density values increases as the excitation energy 
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increases. Caloric curves, i.e. temperature as a function of excitation energy (top panel of 
Fig. [I]), show a monotonic rising behavior. Within statistical errors, a small dependence on 
m, is observed. 
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FIG. 1: (Color online) Temperatures and densities are plotted as a function of the QP excitation 
energy per nucleon for 54 < Aqp < 64. Both quantities are extracted from the quantum method 
based on momentum quadrupole and multiplicity fluctuations. Protons are used as the probe 
particle. Statistical errors are indicated by the bars. 



The density p is plotted as a function of the temperature T in the top panel of Fig. [2] 
for the four different source asymmetries. Universality of a second-order phase transition 
dictates that the different curves in Fig. [2] should collapse to one when the axes are divided 
by their corresponding critical values for each m s . In turn, we can determine the values of 
p c and T c for each m s in order to obtain one universal curve. In the bottom part of Fig. [2] 
we demonstrate that indeed the data displays universality. 
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FIG. 2: (Color online) Top panel: Density plotted as a function of temperature. Bottom panel: The 
reduced density-temperature phase diagram. The dashed line is drawn to guide the eye. Statistical 
errors are indicated by the bars. 



The critical exponent /3 is determined from the relation 



P (^ T 



Pc 



(4) 



This parameter defines the universality class of the system and, therefore, systems with 
similar /3 values have similar underlying physics. Figure |3] shows an excellent fit to the 



data points that results in a slope j3 = 0.35 



expected for the liquid-gas universality class 



IT. 



0.01. This value is in the range of that 



291 ] . The fact that we are getting a value of 



(3 consistent with a liquid-gas phase transition supports our strategy for calculating densities 
and temperatures from the Fermi gas assumption. 
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FIG. 3: (Color online) The extraction of the critical exponent (3, Eq. HI See text for details. 



The pressure of the system is calculated by making use of the grand partition function 
from Fisher's droplet model which is based on the simple idea that a real gas of interacting 
particles can be considered as an ideal gas of clusters (fragments) of various sizes in chemical 



equilibrium |18| . The two-body interaction is assumed to be exhausted in the formation of 
clusters and the resulting clusters behave ideally. The fragment size fluctuations are strongly 
enhanced for systems undergoing a phase transition. The equation of state is expressed in 
terms of the zeroth (Mo) and first (Mi) moments of the fragment (cluster) size distribution 
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TABLE I: Critical values and thermodynamic quantities for the four m s bins. 



m s 


T c (MeV) 


Pc (fm" 3 ) 


P c (MeV/fm 3 ) 


Pc/p c T c 


AH (MeV) 


0.065 


12.12 ± 0.39 


0.070 ± 0.006 


0.211 ± 0.002 


0.25 ± 0.02 


31.50 ± 1.01 


0.115 


12.51 ± 0.35 


0.066 ± 0.005 


0.209 ± 0.001 


0.25 ± 0.02 


32.53 ± 0.90 


0.165 


13.11 ± 0.30 


0.064 ± 0.004 


0.232 ± 0.001 


0.27 ± 0.02 


31.46 ± 0.71 


0.215 


13.39 ± 0.21 


0.061 ± 0.002 


0.258 ± 0.002 


0.31 ± 0.01 


32.13 ± 0.50 



by 

P = T P ^, (5) 

where P is the pressure. The temperature and density have been determined from the proton 
quantum fluctuations and can be used in Eq. [51 The fc-th moment is defined as 

M k = £ A k Y(A) , (6) 

■ro a.r 

where Y(A) is the multiplicity of the fragment A and the largest fragment A max which 
represents the liquid phase is excluded in the summation. The calculated pressure P is 
normalized to its critical value P c and plotted versus T c /T, the inverse of the reduced tem- 
perature in Fig. HI Even though data points from the lowest two excitation energy bins 
deviate from the fitted curve, it is interesting to notice that an excellent fit was obtained in 
the region close to the critical values. The critical pressure is determined using the equation 
.„ PIP, = A-BTJT, with A and B as constants, derived by Guggenheim from the principle 
of corresponding states [30(. Taking A to be nearly equal to B, as it is for van der Waals 
systems, one obtains P/P c = exp[AH /T c (l — T c /T)] which is the Clausius-Clapeyron equa- 
tion that describes several fluids for T < T c . The quantity AH is the enthalpy of emission 
of a fragment from the liquid. 

The experimental critical parameters (T c , p c and P c ) that give the location of the critical 
point along with the critical compressibility factor for each m s bin are listed in Table [H 
The critical temperatures T c are observed to increase when increasing m s . This trend is 



consistent with the results reported in Ref. 



25] . On the other hand, the critical densities 



p c decrease with increasing m s while the critical pressures P c generally increase with m s . 
The critical compressibility factor values (P c /p c T c ), which quantify the deviation from ideal 
gas, in turn increase when increasing m s . These values are very close to those of real gases 
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FIG. 4: (Color online) The reduced pressure as a function of the inverse of the reduced temperature 



301 ] according to the principle of 



(symbols with statistical error bars) as derived by Guggenheim 
corresponding states. The solid line shows a fit to the Clausius-Clapeyron equation where the value 
of 2.5 is the average of the AH/T C values extracted from different m s bins. 

0, 0I- The values of AH are also reported in Tableland are observed to weakly depend 
on m s . 

By means of a new quantum method and the grand partition function of Fisher's droplet 
model, the temperature, the density and the pressure of the selected fragmenting sources 
have been calculated and the phase diagrams have been constructed. These parameters and 
the corresponding critical values have shown a dependence on the source neutron/proton 
concentration. A complete experimental location of the critical point is given. The critical 
exponent (3 and the critical compressibility factor have been extracted and found to belong 
to the liquid-gas universality class. Strong evidence for a signature of a liquid-gas phase 
transition in two-component systems has been found. These results provide a means to 
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establish the proton-fraction dependence of the EOS in systems with large neutron excess 
such as neutron stars. 
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